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Abstract. We review some recent results of the theory of Lie systems in order to apply 
such results to study Ermakov systems. The fundamental properties of Ermakov systems, 
i.e. their superposition rules, the Lewis-Ermakov invariants, etc., are found from this new 
perspective. We also obtain new results, such as a new superposition rule for the Pinney 
equation in terms of three solutions of a related Riccati equation. 
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1 Introduction 

Nonlinear equations have been of an increasing interest in Physics during the last thirty years 
and one of the simplest examples is the today called Milne-Pinney equation 

x = -u 2 (t)x + -=, (1) 
x° 

where k a is real constant with values depending on the field in which the equation is to 
be applied. Such equation was introduced by an Ukrainian mathematician of the nineteenth 
century, V.P. Ermakov, as a way of looking for a first integral for the time-dependent harmonic 
oscillator pQ. He was using some of Lie's ideas for dealing with ordinary differential equations 
with tools of classical geometry. Lie had obtained a characterisation of non-autonomous systems 
of first-order differential equations 

— = Y*(t,x), i = l,...,n, (2) 

admitting a superposition rule [2]. Such a problem has being receiving these last years very 
much attention because of its very important applications in physics and mathematics, see [31 
HI [51 El El El HUl HI H21 HSl HH HSl [El [IT]. This approach has recently been revisited from a 
more geometric perspective in [18] where the role of the superposition function is played by an 
appropriate connection. This new approach allows us to consider a superposition of solutions 
of a given system in order to obtain solutions of another system as it was done in [19]. Our 
aim in this paper is to show how such a superposition rule may be understood from a geometric 
viewpoint in some interesting cases, as the Milne-Pinney equation ([1]) [20\ I21j . the Ermakov 
system and its generalisations. 

Recall that Ermakov systems are systems of second-order differential equations composed 
by the Milne-Pinney differential equation ([!]) together with the corresponding time-dependent 
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harmonic oscillator. They have received lot of attention since its introduction in Ermakov's 
paper pQ, where many of its important properties have been exhibited. For instance, he found the 
time-dependent invariant called Ermakov-Lewis invariant. In a short paper Pinney [21] was able 
to prove the fundamental result that the general solution of the Milne-Pinney equation can be 
described in terms of an arbitrary pair of independent particular solutions of the corresponding 
time-dependent harmonic oscillator and two constants. 

Ermakov systems have also been broadly studied in Physics since its introduction until now. 
They appear in the study of the Bose-Einstein condensates and cosmological models [22 ] [23 | [24] . 
and in the solution of time-dependent harmonic or anharmonic oscillators |25[ [26] [27] [28] [29] . 
Much work has also been devoted to use Hamiltonian or Lagrangian structures in the study of 
such a system, see e.g. [30} 131] . and many generalisations or new insights from the mathematical 
point of view can be found in[32l[Ml[3l[35l[Ml[371[Ml[39l[10]. 

Our aim in this paper is to review the mathematical insight developed in [6] in order to 
study the properties of Ermakov systems. With this objective, we review in Section [2] the 
concept of superposition rule, a map describing the general solution of a certain kind of system 
of differential equations called Lie systems, in terms of a finite set of particular solutions of it. 
The concept is generalised in Section[3]to express the general solution of a system of second-order 
differential equations (hereafter shortened as SODE) in terms of arbitrary particular solutions of 
a new Lie system and some constants. We prove that the theory developed by Lie can easily be 
adapted to deal with many examples of such SODE systems, and in this way the superposition 
rules which are generally used for first-order differential equations can also be used for dealing 
with Ermakov systems. We find room in this way for implicit nonlinear superposition rules 
in the terminology of |36[ 141] . The Ermakov-Lewis invariants [29] appear in a natural way as 
functions defining the foliation associated to the superposition rule. Some simple examples are 
used to illustrate the theory. Section 0] is devoted to describe the possibility of a superposition 
of solutions of a system to give a solution of another related one, and we show how Pinney's 
results [21] fit in this framework, i.e. we recover the known superposition rule in which the 
solutions for the Pinney equation are expressed in terms of two solutions of the corresponding 
i-dependent harmonic oscillator and as a straightforward application of our method we show 
that the general solution of the Pinney equation can also be written as a superposition of three 
particular solutions of the Riccati equation. 

2 Systems of differential equations 
admitting a superposition rule 

A superposition rule for solutions of ([2]) is determined by a function $ : R n ( m+1 ) — > M. n , 

x = ^{x^),...,x^ m y,ki,...,k n ), x (a) e R n , a = l,...,m, 
such that the general solution can be written, for sufficiently small t, as 

x(t) = $(x(i)(*),...,X( m )(i); 

where {x( a \(t) \ a = 1, ... , m} is a fundamental set of particular solutions of the system ([2]) 
and k = (hi, . . . , k n ) is a set of n arbitrary constants, whose values fix each particular solution. 
A superposition rule provides us a method for finding the general solution of ([2]) in terms of 
a set of particular solutions. 

According to the Implicit Function Theorem, the function $ : ]R n ( m+1 ) — > W 1 can be, at least 
locally around generic points, inverted, so we can write 

k = *(S(0),---)Z(m)) 
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for a certain function ^ : R n ( m+1 ) — > R n . Hereafter in order to handle a shorter notation 
we started writing xr \ instead of x. The level sets of the superposition function ^> define 
a foliation that is invariant under permutations of the (m + 1) variables. The functions 3> and 
are related by: 

k = . . ,X( m y,ki,. . .,k n ),x {1) ,. . . ,x (m) ). (3) 

The fundamental property of the superposition function ^ is that as 

* = ^(x( )(t),x {1) (t),...,x {m) (t)), 

the function *(z(o), • • • , X( m )) is constant along any (m + l)-tuple of solutions of the system ([2]). 
This implies that if the 'diagonal prolongation' Y(t, X( ), • • • , X(m)) °f the t-dependent vector 
field Y(t,x) = Y*(t,x)d/dx t is defined by 

m 

Y(t,x {0) ,... ,x (m) ) = ^Y a {t,x {a) ), t G R, 



a=0 



where 



n d 
Y a (t,x {a) ) = J2Y t (t,x (a) " 



dx) , 

i=l w 

then Y(t, X( ), • • • , X( m )) is a t-dependent vector field on R n ( m+1 ) tangent to the level sets of 'J/, 
i.e. the components *$> l are constants of motion. The level sets of ^ corresponding to regular 
values define a n-codimensional foliation T on an open dense subset U C R n ( m+1 ) and the family 
{Y(t, •), t G R} of vector fields in R n ( m+1 ) consists of vector fields tangent to the leaves of this 
foliation. 

Remark that, as pointed out in [TS], for any k = (ki, ■ ■ ■ , k n ) G R n and each (xm, . . . , Xi m \) G 
R nm there is one point (x( ),X(i), • • • >£( m )) on the level set of the foliation !F, namely, such 
that . . . , xi m )\ k), X(i\, . . . , xi m )) ^ ( c f- ©); then, the projection onto the last m 

factors 

pr : (x (0) ,x (1) , . • .,x (m) ) G R<™ +1 ) » (x (1) , . . . ,x (m) ) G R nm 

induces diffeomorphisms on the leaves J~k of F. Such a foliation gives us the superposition 
principle without referring to the function ^f: if we fix the point xm\(0) and m solutions 
. . . , X( m )(t), then X(o)(i) is the unique point in R n such that (xr \(t),X(^(t), . . . , X( m )(i)) 
belongs to the same leaf of T as (x( O )(0),xn)(0), . . . , X( m )(0)). This means that only F really 
matters for the superposition rule. 

Lie's main result [2] can be expressed as follows: 

Lie's Theorem. The system ^ on a differentiable manifold N admits a superposition rule if 
and only if the t-dependent vector field Y(t, x) can be locally written in the form 

r 

Y(t,x) = J2 b a(t)X a (x), (4) 

a=l 



where the vector fields X a , with a = 1, . . . , r, close on a r -dimensional real Lie algebra, i.e. there 



exist r 3 real numbers c„« 7 such that 



[X a ,Xp] = ^2c a/3 ~ ( X 1 , Va,f3 = 1, . . . ,r. (5) 

7=1 

In this case we say that ^ is a Lie system. 
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The number m of solutions involved in the superposition rule for the Lie system defined 
by (UD with generic b a (t) is the minimal k such that the diagonal prolongations of X\, . . . ,X r 
to N k , X±, . . . , X r , are linearly independent at (generically) each point: for k = m the only real 
numbers solution of the linear system 

r 

^2c a X a (x w ,...,x {k) ) = 0, 

a=l 

at a generic point (xm, . . . , X(ia) is the trivial solution c a = 0, with a = 1, . . . , r, and there are 
nontrivial solutions for k < m. Then, the superposition function ^ is made up of t-independent 
constants of motion for the prolonged i-dependent vector field Y. We shall call them first- 
integrals. 

A possible generalisation consists on considering foliations of codimension different from n, 
or even choosing different Lie systems with the same associated Lie algebra for defining the 
prolonged vector field. These facts will be illustrated with several examples in the following 
sections. 



3 SODE Lie systems 

A system of second-order differential equations 

d 2 x % , / dx\ 

F % [ tx, — , i = l n 



dt 2 V dtj' 

can be studied through the system of first-order differential equations 

x % = v\ 

v l = F l (t,x,v), i = l,...,n, 
with an associated t-dependent vector field 

We call SODE Lie systems those SODE for which X is a Lie system, i.e. X can be written as 
a linear combination with t-dependent coefficients of vector fields closing a finite-dimensional 
real Lie algebra. Even if we are focusing here our attention on second-order differential equations, 
the approach can be straightforwardly applied for r-order differential equations. Some particular 
examples of second-order Lie systems are given in next subsections. 



3.1 Milne Pinney equation 

The Milne-Pinney equation ([T]) is a second-order nonlinear differential equation [20l [21] that 
describes the time-evolution of an isotonic oscillator [121133] (also called pseudo-oscillator), i.e. 
a harmonic oscillator with inverse quadratic potential [32]. This nonlinear oscillator shares 
with the harmonic one the property of having a period independent of the energy [H], i.e. it is 
a isochronous system, and its quantum version has an equispaced spectrum [4"5| 135]. 

Note that if x(t) is solution of ([1]) then x'(t) = —x(t) is solution too and we can restrict 
ourselves to consider solutions in the half-line R + = {x 6 R | x > 0}. Thus, we are interested in 
the solutions of ([!]) that are curves in R + . 
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We can relate the Milne-Pinney equation with a system of first-order differential equations 
defined in TR + by introducing a new auxiliary variable v = x. Such system is given by 



k 

u = -w 2 (t)a; + — , (6) 
x 6 

where x G R+ and (x, v) G T X R + . By the sake of simplicity we assume from now on that k > 
but similar results can be obtained for k < 0. Now, the associated t-dependent vector field 
over TR + is 

d f 2,, k \ 9 
X = v— +[ -uj\t)x + -»)-£-. 
OX \ X 6 J Ov 

This is a Lie system because X can be written as 
X = L 2 -uo 2 {t)L 1 , 

where the vector fields L\ and L2 are given by 

d d k d 

OV OX X 6 ov 

which are such that 

[L\, L2] = 2L3, [L3, L2] = — L2, [Ls,Li]=Li 

with 

1/5 d 

L3 = — x— v— 

2 \ ox ov 

i.e. they span a 3-dimensional real Lie algebra g isomorphic to s((2,R). Another different 
relation of the Lie algebra s[(2,M) with this system has been pointed out but in the context 
of Lie symmetries |471 148] , On the contrary here s[(2, M) is the Lie algebra of vector fields 
of the Lie system and the search for the general solution of (0) can be reduced to finding 
the solution starting from the neutral element of the equation in the group SL(2, R) given by 
g g~ l = u; 2 (t)ai — a,2 (see [5j[6]), where {ai,a2,aa} is the following basis of st(2,R): 

/ \ / -1 \ 1 / -1 \ 

ai = { -1 ) ' a2 = { ) ' &3 = 2 { 1 ) ' (7) 

which satisfy the same commutation relations as the vector fields {Lj \ j = 1, 2, 3}. Actually, it 
is possible to show that Lj is the fundamental vector field corresponding to &j with respect to 
the action $ : (A, (x,v)) G SX(2,R) x TR + -» (x,v) G TR + given by: 



{(3v + ax){5v + 7a;) + k{8(5/x 2 ) 




(Sv + 7x) 2 + k(5/x) 2 



where k is ±1 or depending on the initial point (x,v) and the element of the group SL(2,M) 
that acts on it. In order to obtain an explicit expression for k in terms of A and (x,v) we can 
use the following decomposition for any element of the group SL(2, R) 



A = exp(— a2&2) exp(— aiai) exp(— c^as) 



i 


OL2 \ 


1 


! ( 




1 J 


V oil 





e « 3 /2 o 
e- Q3 / 2 
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from where we obtain that a\ = j6 and «2 = P/S. As we know that <3?(exp(— 0^3), (x,v)) is 
the integral curve of the vector field L 3 starting from the point (x,v) parametrised by 03, it is 
straightforward to check that 

(xi,vi) = $>(exp(-a 3 a 3 ), (x,v)) = (exp(a 3 /2)x, exp(-a 3 /2)v), 

and in a similar way 

(x 2 ,v 2 ) = $(exp(-aiai), (a;i,ui)) = (xi,a±xi +v{). 

Finally, we want to obtain (x, v) = $(exp(— 02^2), (x 2 , v 2 )), and taking into account that the 
integral curves of L2 satisfy that 

x 3 dv dx 

— = - = da 2 , (8) 

it turns out that when k > we have v 2 + k/x 2 = v 2 + k/x 2 = A with A > 0. Thus, using this 
fact and ([8]) we obtain 

k l / 2 dv 

dao, 



{\-v 2 fl 2 

and integrating v between V2 and v, 

V A v 2 1 , . , , N 

«2TTrt + 7^ = TT77 (^A + v 2 \x 2 \) ■ 



(A — v 2 ) 1 / 2 "k 1 / 2 (A -v\yi 2 k 1 / 2 
As k = sign[v], we see that k is given by 



sign[a 2 A + v 2 \x 2 \] = sign 



{3 . r , kSP \x\ . . 

T (x7 + vo) H 7T + -r [vd + x^y) 

d x z 



Non-trivial first-integrals independent of u(t) for the equation ([T]) do not exist, i.e. there is not 
any i-independent constant of motion I : U C TM + — > R such that XI = for any function u)(t). 
This is equivalent to dl(Lj) = for Lj with j = 1,2,3. Thus, the integrals of motion we are 
looking for verify that dl vanishes on the involutive distribution V(x) ~ (L\(x), L 2 (x), L^(x)) 
generated by the fundamental vector fields Lj. In almost any point we have that V{x) = T X TIR + , 
and as dl = on dense subsets in TR+, the only possibility is dl = 0, therefore / is a constant, 
and the integral of motion is trivial. 

3.2 Generalised Ermakov system 

Consider a possible generalisations of the Ermakov system given by 



x = -u; 2 (t)x + -^f(y/x), 
x A 

1 



X 

2 



y = -u\t)y + —g{y/x). (9) 

This system of differential equations has been broadly studied in [271 ESS EHl EH HU HH1 [50]. In 
this section we analyse how can we make use of the theory of Lie systems in order to obtain the 
known first-integral of motion for this system of differential equations. 

Note that if (x(t),y(t)) is a solution, then (— x(t), —y(t)) is also a solution for Q, and therefore 
we can focus our attention on the solutions of generalised Ermakov system which are curves in 
the manifold M+ x R. 
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The generalised Ermakov system ([9]) can be written as a system of first-order differential 
equations with a double number of degrees of freedom by introducing two new variables v x 
and v y , 



x = v a 



-cj 2 (t)x + ^f(y/x), 



V = Vy, 



X 3 ' 

v y = -uJ 2 (t)y + ^g(y/x), (10) 

y 6 

and its solutions can be studied from the integral curves of the t-dependent vector field X in 
T(R+ x R) given by 

x = v * Tx + | + {~^ t)x + lk + ("" 2(0y + ¥ 9{y/x) ) ~k> 

which can be written as a linear combination 
X = N 2 -uj 2 (t)N 1 , 

where N± and N 2 are the vector fields 

d d d d 1 d 1 d 

N\ = x- h y— , iV 2 = v x — + v y — + -»f(y/x)- h -zg{y x) — . 

ov x ov y ox oy x A ov x y 6 av y 

Note that these vector fields generate a 3-dimensional real Lie algebra with a third generator 

AT -- ( — —- — - — 

3 2 \ dx ^ dy Vx dv x Vy ' dv y 

In fact, as 

[N lt N 2 ] = 2N 3 , [N 3 , Nt] = iVi, [N 2 , N 3 ] = N 2 , 

they generate a Lie algebra isomorphic to st(2,R). Therefore the generalised Ermakov system 
is a SODE Lie system with associated Lie algebra sl(2, R). 

The integrable distribution associated to the Lie system (|10[) is of dimension three, while the 
manifold T(R + xR) is 4-dimensional, and then there exists a first-integral of motion, F : R 4 — > R, 
for any to 2 (t). This first-integral F satisfies N { F = for i = 1,2,3, but as [iVi,^] = 2N 3 it 
is enough to impose N%F = N 2 F = 0. Then, if N\F = 0, by the method of characteristics we 
can conclude that there exists a function F : R 3 — ► R such that F(x, y, v x ,v y ) = F(x, y, £) with 
^ = xv y — yv x ). The condition N 2 F = reads now 

dF dF ( y r . , , i , . A flf 

We can therefore consider the associated system of the characteristics: 

dx dy d£ 

v x v y -% s f(y/x) + ^g(y/x)' 

and using that 

—ydx + xdy dx dy 

i V X Vy ' 
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we arrive to 

—ydx + xdy d^ 



i.e. 



£ ~ X s fix ) + §s9ix) 

and integrating we obtain the following first-integral 



M?) +< «(? 



dC = C. 



This first-integral allows us to determine, by means of quadratures, a solution of a subsystem 
in terms of a solution of the other equation. 

3.3 The harmonic oscillator with time-dependent frequency 

Our aim in this section is to exhibit with a simple example that we can define new Lie systems 
of differential equations admitting first-integrals of motion which do not depend on the time- 
dependent coefficients by putting together some copies of the same Lie system in the way 
indicated in Section [2j These integrals of motion allow us to obtain relations between different 
solutions of the initial Lie systems and can be used to construct a superposition rule if enough 
copies have been added. 

The equation of motion of a 1-dimensional harmonic oscillator with time-dependent frequency 
is x = —oj 2 (t)x, which is a second-order differential equation whose solutions are curves in R. 
We can alternatively consider the following system of first-order differential equations 

x = v, 

v = -uj 2 (t)x, (11) 
whose solutions are the integral curves in TR of the t-dependent vector field 
d 9,,n d 

Now X is a linear combination X = X 2 — u> 2 (t)Xi, where X\ and X2 are vector fields in TR: 

Y 9 V 9 

Xi = x—, X 2 = v—, 
ov ox 

such that 

[X 1 ,X 2 ] = 2X 3 , [X 3 ,X 1 ]=X 1 , [X 2 ,X 3 }=X 2 , (12) 
where X3 is the vector field given by 

2 \ dx dv 

As a consequence of Lie's Theorem we see that X defines a Lie system with associated Lie 
algebra sl(2,R). Actually, the vector fields {X a \ a = 1,2,3} are fundamental vector fields 
corresponding to the usual linear action of SL(2, R) on R 2 ~ TR given by 



x = ax + f3v, v = 7X + Sv if A 



a j3 
7 5 
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Remark that the elements of the basis of st(2,R) given by ([7]) close the same commutation 
relations as the vector fields {X a \ a = 1, 2, 3} which turn out to be the associated fundamental 
vector fields for the elements & a . 

Non-trivial first-integrals independent of uj(t) for the system (|lip do not exist, because the 
fundamental vector fields generate the involutive distribution V — (Xi, X2, X 3 ) whose values 
in each point of TM generate the corresponding tangent space in almost any point and then 
as dl vanishes on V, therefore dl vanishes in a generic point. Thus / is a trivial first integral of 
motion. 

As it was said before, non-trivial integrals of motion are needed to develop the procedures 
of integration indicated in [6], i.e. we should add copies of the same differential equation until 
the linear subspace of the involutive distribution of the total system in a generic point has 
a dimension smaller than the dimension of the manifold. In this way there is room for a non- 
constant function / such that dl vanishes on the subspace of the distribution associated with 
the extended system of differential equations. 

Thus, we first consider a system with two copies of the same differential equation, 

'±\ = —uj 2 (t)xi, 

'±2 = -uj 2 (t)x2, (13) 

which corresponds to a 2-dimensional isotropic harmonic oscillator with a time-dependent fre- 
quency and is associated with the following system of first-order differential equations 

xi = Vl, 

v\ = —uj 2 (t)xi, 

£2 = v 2 , 

V 2 = -L0 2 (t)x 2 , (14) 

whose solutions are the integral curves of the t-dependent vector field 

d d 2 . x d 2 , . d 

X = Vl- h V 2 -z U) {t)Xi U) (t)X2 ^— , 

OXi OX2 OVi OV2 

which is a linear combination, X = X2 — u> 2 (t)Xi, with Xi and X 2 being the vector fields 

d d d d 

Xi = xi- hx 2 ^— , X 2 = vi- \~v 2 -^ — , 

OVi OV2 0x1 0x2 

which are such that 

[Xi,X 2 ] = 2X 3 , [X 3 ,X 1 ]=X 1 , [X 2 ,X 3 }=X 2 , (15) 
where the vector field X 3 is defined by 

1 / d d d d 

X 3 = o ^^2^ Vl- V 2 ^— 

2 V oxi 0x2 ovi 0V2 



The system of differential equations (|14p is therefore a Lie system, i.e. the t-dependent vector 
field can be written as a linear combination with time-dependent coefficients of vector fields X a 
closing a finite dimensional Lie algebra. Such vector fields are diagonal prolongations of the 
fundamental vector fields of the initial system of differential equation and satisfy the same 
commutation relations (fT5|) as (fT2|) . The subspace of the distribution associated for this Lie 
system has a rank lower or equal to the dimension of the Lie algebra. The Lie algebra associated 
with the t-dependent vector field X is also sl(2, R). 
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Note that the system (jlip did not admit non-trivial first-integral of motion independent 
of uj(t), but the system (fP5]) does admit a first-integral because the subspace of the distribution 
generated by the fundamental vector fields has rank three in almost any point and the dimension 
of the total manifold is four. Thus, there exists a function F such that dF is vanishes on the distri- 
bution generated by the X a 's. The function F(x±, X2,vi,V2) is such that X\F = iff there exists 
a function F{x\, X2, £) with £ = x\V2 — X2V1, such that F(xi, X2, v±, V2) = F(xi,x 2 , £), an d then 
the second condition X2F = implies that F can only depend on £, i.e. F(xi,x 2 ,vi,v 2 ) = F(£). 
From the commutation relation 2X3 = [Xx,^], we see that the conditions X\F = X2F = 
imply X3F = 0, and therefore any first integral is a function of F(x±, X2, v%, V2) = x\V2 — X2V1, 
which physically corresponds to the angular momentum. This first integral can be seen from 
the mathematical viewpoint as a partial superposition rule [6]. Actually, if x\{t) is a solution 
of the first equation, then we obtain for each real number k the first-order differential equation 
for the variable X2 

. . dx2 , . / n 

= k + xi(t)x 2 , 

from where X2 can be found to be given by 



x 2 {t) = k'xi(t) + kxi{t) 



where k' is a new integration constant. 

In order to look for a superposition rule we should consider a system of some copies of (jlip 
with at least as many integrals of motion as the dimension of the initial manifold. Then it may be 
possible to obtain the variables of the initial manifold explicitly in terms of the other variables. 
Following the development of the Section [2] a way to obtain the number of necessary particular 
solutions to obtain a superposition rule is to consider a set of copies of the initial system and 
check out if the prolongations of the vector fields X\, X 2 and X3 are linearly independent in 
a generic point. 

In the case of two copies of the t-dependent harmonic oscillator it is possible to show that in 
a generic point if \±X± + X 2 X 2 + A3X3 vanishes then Ai = A2 = A3 = 0, therefore m = 2 and 
consequently there is a superposition rule involving two particular solutions. In order to obtain 
the superposition rule we need to consider three copies of the ^-dependent harmonic oscillator 
and study the system of first-order differential equations 

±1 = Vi, 

vi = -uj 2 (t)xi,x 2 = v 2 , 
v 2 = -uj 2 (t)x2, 
x = v, 

V = —L0 2 (t)x, 

whose solutions are the integral curves of the t-dependent vector field 

d d d 2, \ d ii \ 9 ^ 1 \ d 

X = V 1T - \-V 2 -^ hV- UJ {t)Xl- LO {t)X2T ) W {t)x — , 

OX 1 OX2 OX OV\ ov 2 ov 

which is a linear combination, X = X 2 — L0 2 (t)Xi, with X\ and X2 being the vector fields 

d d d d d d 

Xi=Xi- VX 21 . X 2 =V 1 - h«27 ^ V ^~i 

OV\ OV2 OV OXl OX2 OX 

which are such that 



[Xl,X2] — 2X3, [X3,Xi] — Xi, 



[X2, X 3 ] — X 2 , 
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where the vector field X3 is defined by 

X - - ( — — — — — —\ 

3 2 \ 1 dxi X2 8x2 X dx Vl dvi 2 dv2 V dv J ' 

The first-integrals F are as in the last case the solutions of X\F = X2F = 0, because 
2X3 = [X 1 ,X 2 \. The condition X\F = says that there exists a function F : R 5 -> M 2 
such that F(x±,X2,x, v±,V2,v) = F(x\,X2,x, £1, £2) with £i(xi, X2, x, v\,V2,v) = xv\ — x\v and 
£2(^1, %2, vi,V2,v) = XV2 — X2V, and when written in terms of F the condition X2F = implies 
that F is an arbitrary function of £1 and £2 (of course, £ = X\V2 — X2V1 is also a first-integral). 
They produce a superposition rule, because from 

XV2 — X2V = k\, X\V — V\X = k2 

we obtain the expected superposition rule for two solutions: 

h 

X = C±Xl + C 2 X2, V = ClVl + C 2 V2, Ci = — : K = X\V2 — X2V\. 

k 

4 Construction of mixed superposition rules 

In the preceding sections we have considered SODE Lie systems together with some copies of 
the same system. In this section we show a way to deal with sets of maybe different systems of 
differential equations sharing the same associated Lie algebra as Lie systems. These sets satisfy 
certain conditions to be explained later on and can be used to obtain mixed superposition rules. 
Suppose that a Lie system is given 

where the X a with a = 1, . . . , r, close a finite-dimensional Lie algebra as given by (|5|). We want to 
obtain a mixed superposition rule for this Lie system, i.e. a function $ : iV 1 x • • • x N m x ]R n — > N 
in such a way that any solution integral curve x(t) of X(t, •) be given by 

x(t) = $(x(i)(i) } ... ,x m (t),kx,... ,k n ), 

where x^(t), . . . , #( m )(t), are integral curves of the £-dependent vector fields on the correspon- 
ding manifolds N 3 with j = 1, . . . , m and given by 

r 

X^(t,.) = Y,b a (t)X 1 a (-), 
a=l 



xM(t,.) = 5> a (t)JTO 

a=l 

and in such a way that Xa with j = 1, . . . , m, and a = 1, . . . , r, close the same commutation 
relations as the vector fields X a , i.e. 

[X 3 a ,X 3 p] = c^X^ j = l,...,m, and a, /3 = 1, . . . r. 

The search for a set of vector fields closing a given Lie algebra in a certain manifold amounts 
to look for different actions of the corresponding Lie group. 
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Then, over N we obtain that 

[X a ,X 3 p]=0, [Xi,Xj$)=0, j,k = l,...,m, with j ^ k, and a,0 = l,...r, 

and because of these commutation relations we can define in TV the vector fields 

m 

Ya = X a + ^2 X„ 

i=l 

which satisfy the same commutation relations as the vector fields X a in ([5]). Thus, the system 
of differential equations determining the integral curves of the t-dependent vector field 

r 

a=l 

is a Lie system with the same associated Lie algebra as the initial system of differential equations. 
Nevertheless, the dimension of the manifold is this time larger than the dimension of iV because 
we have put together several manifolds. The distribution V associated to this Lie system is 
given by V{x) = (Y"i(x), . . . ,Y r (x)) and is involutive. The dimension of the subspace in a point 
is lower or equal to r and therefore when the dimension of the manifold N is larger than r there 
are integrals of motion for any value of the time-dependent coefficients which may be used to 
obtain superposition rules. 

In the next section we give some examples of the application of this procedure to the Ermakov 
system in order to both recover previously known properties and to find also new results. 



4.1 Ermakov system 

Consider the system of ordinary first-order differential equations 



X — V x , 

,2 



v x = to (t)x, 
y = Vy, 



1 



made up by a 1-dimensional harmonic oscillator and the Milne-Pinney equation with k = 1, 
which has the above-mentioned conditions. Its solutions are the integral curves of the i-depen- 
dent vector field 

d d o, s d f o,., 1\ 3 



which is a linear combination with time-dependent coefficients, X = X2 — uJ 2 {t)X\, of 
d d „ d did 



x i = x x h ^« — ' X 2 = v x — + v y —+ . 

dv x dv y ox dy y- 3 dv y 

This system closes on a sl(2, M) Lie algebra with X% given by 

1/(9 d d d 



^ 3 2 \ dx ^ dy Vx dv x Vy dv y 

The generators of this Lie system span a distribution of dimension two and there is no 
first-integral of the motion for such subsystem. By adding the other sl(2,R) linear Lie system 
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appearing in the Ermakov system, the harmonic oscillator with the same time-dependent angular 
frequency, as the rank of the space of the distribution in the 4-dimensional space is three, there is 
an integral of motion. The first-integral can be obtained from X\F = X 2 F = 0. But X\F = 
means that there exists a function F : K 3 — > M such that F(x,y,v x ,v y ) = F(x,y,£), with 
£ = xv y — yv x , and then X2F = is written 

dF dF x dF 
x dx v dy y 3 d£ 

and we obtain the associated system of characteristics 

xdy - ydx _ y 3 d£ d(x/y) yd^ _ 
£ x £ x 

from where the following first-integral is found [29]: 



X \ .0 / oc 



2 



ip(x,y,v x ,v y ) = ( - J +£ = (-1 +(xv y -yv x ), 
which is the well-known Ermakov-Lewis invariant [27^ [28| HO] . 



4.2 The Pinney equation revisited 

We mentioned before the possibility of obtaining solutions of a given system from particular 
solutions of another related system. We next study a particular example which allows us to 
recover the results obtained by Pinney long time ago [21]. Consider the system of first-order 
differential equations: 

X — V X i 

y = v y , 

z = v z , 

k 

V x = -L0 2 {t)x + -3, 

Vy = -u?(t)y, 

which corresponds to the vector field 

d d d k d 2/ \ f ® ® ® \ 
x dx Vy dy Vz dz x 3 dv x \ dv x ^ dv y Z dv z J 

The t-dependent vector field X can be expressed as X = N2 — u> 2 (t)Ni where Ni and N2 are 

d d d d d did 

m=y- Vx- Vz-—, N 2 = v x — +v y — + v z — + , 

dv y dv x dv z dx dy dz x 6 dv x 

These vector fields generate a 3-dimensional real Lie algebra with the vector field -/V3 given by 

N --( — — —- — - — - 8 \ 

3 2 y dx ^ dy 2 dz x dv x y ' dv y Vz dv z ) ' 

In fact, they generate a Lie algebra isomorphic to st(2,R) because 



[Nx, N 2 ) = 2N 3 , [N 3 , Ni] = iVi, [N 2 , N 3 ] = N 2 . 
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The dimension of the distribution generated by these fundamental vector fields is three 
and the manifold of the Lie system is of dimension six, then there are three time-independent 
integrals of motion which turn out to be the Ermakov invariant I\ of the subsystem involving 
variables x and y, the Ermakov invariant I2 of the subsystem involving variables x and z, and 
the Wronskian W of the subsystem involving variables y and z. They define a foliation with 
3-dimensional leaves. This foliation can be used for obtaining a superposition rule. 

The Ermakov invariants read as 

h = 2 ( (y v x - xv y ) 2 + k j , h = g ( ( xv z ~ zv x) 2 + k (^j 

where I\ and I 2 are non-negative constants and the Wronskian W is: 

W = yv z - zv y . 

We can obtain an explicit expression of x in terms of y, z and the three first integrals I\, I2, W: 

x = ^Q- [hy 2 + hz 2 ± \f±hh - kW 2 y Z y /2 . 

Remark that here W is a constant fixed by the two independent particular solutions of the 
time-dependent harmonic oscillator x\{t) and X2(t) and only I\ and I2 play the role of constants 
in this superposition rule for the Milne-Pinney equation. This is not a surprising fact because 
the Milne-Pinney equation is a second-order differential equation. Note also that the values 
of I\ and I2 are non- negative constants but should be chosen such that x(0) be real. 

This can be interpreted, as pointed out by Pinney [21], as saying that there is a superposition 
rule allowing us to express the general solution of the Milne-Pinney equation in terms of two 
independent solutions of the corresponding harmonic oscillator with the same time-dependent 
angular frequency. 

4.3 A new superposition rule for the Pinney equation 

A new mixed superposition rule involving solutions of a Riccati for obtaining the general so- 
lution of Pinney equations is obtained in this section as a straightforward application of our 
development. 

The t-dependent Riccati equation 

^ = b (t) + b 1 (t)x + b2(t)x 2 . (16) 

has been studied in [5 1 1. [52] from the perspective of the theory of Lie systems. We follow 
here a very similar approach. From the geometric viewpoint the Riccati equation (I16j) can 
be considered as a differential equation determining the integral curves of the time-dependent 
vector field 

r = (b (t) + b 1 (t)x + b 2 (t)x 2 )^-, 

which is a linear combination with time-dependent coefficients of the three vector fields 

t 9 t 2 9 r 9 

Li = -z-, L 2 = -x — , L 3 = -x—, 

ox ox ox 

which close on a 3-dimensional real Lie algebra, with defining relations 

[L 1 ,L 2 ] = 2L 3 , [L 3j L 1 ]=L 1 , [L 2 ,L 3 ]=L 2 , (17) 



therefore isomorphic to s[(2, M), because the commutation relations (fTTj) are the same as (fl~2j) . 
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The following particular case of Riccati equation: 



dx 
~dt 



x 



is the equation for the integral curves of the t-dependent vector field L = L2 — u> 2 (t)hi. Con- 
sequently, we can apply the procedure of the Section H] and consider the following system of 
differential equations 



±1 = — oJ 2 (t) 
±2 = —oj 2 (t) 
±3 = -u?(t) 
x = v, 



r 2 
x l5 

T 2 
x 2i 

T 2 



v = —uj 2 (t)x H — ^ 

X 6 

described by a vector field in R 3 x TR + , where (xi,x 2 ,X3) € R 3 , x £ 
According to our general recipe, consider the following vector fields 



and (x, v) £ T x 



d d d d 
ax\ 0x2 0x3 ov 



Ah 















-xi- 



X2- 



X3- 



M 2 = -x\ A 

OX l 

ltd d 

+ ~ 1 x a 

OX ov 



,_d_ 

l dx, 







+ v— + — — 

3 dx-i dx x 3 dv ' 



dx\ 8x2 dx% 2 
that, by construction, verify the same commutation relations as before, i.e. 
[Mi,M 2 ] = 2M 3 , [M 3 , Mi] = Mi, [M 3 , M 2 ] = -M 2 , 

and the full system of differential equations can be understood as the system of differential 
equations for the determination of the integral curves of the t-dependent vector field M = 
M2 — uj 2 (t)Mi. The dimension of the distribution associated to this Lie system is three in 
almost any point and then there exist two integrals of motion. As 2M 3 = [Mi, M 2 ], it is enough 
to find the simultaneously integrals of motion of Mi and M 2 , i.e. a function F : M 5 — ► M such 
that M\F = M 2 F = 0. 

Let us first look for first-integrals independent of x 3 , i.e. we suppose that F depends just 
on xi, x 2 , x and v. Using the method of characteristics, the condition Mi-F = means that 
there is a function F : M 3 — > R such that F(x±, x 2 , x, v) = F(Ii, h, h) with I\, I 2 and I3 given by 

Ii=xi-x 2 , I 2 = X2-v/x, h = x. 

Then, the condition M2F = reads in terms of the variables I\, I 2 , I3 and I± = v as 



h 



h dF I 2 OF dF 

-2— 2— 1 

I 3 dh h dl 3 dh 



+ 



-I( - 2hl 2 



dF 
dh 



I 2 2 + 



dF 
dh 



0. 



Thus, the linear term in h and the other one must vanish independently. The method of 
characteristics applied to the first term means that there exists a function F : R 2 — > R such 
that F(h,h,h) = F{K 1 ,K 2 ) where 



h 
I2 



K 2 



hll 



Finally, taking into account the last result in M 2 i ? = we obtain 



-A 2 - Ax + 



kKA dF 



K 2 
^2 



dKx 



A 2 + 



k 



dF 



A 2 J dK 2 



0. 
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and by means of the method of characteristics 
d Kl _ Kl + K l - k -^f 

Finally, taking into account the last result we obtain the first-integral 



Ci = K 2 + 



k + Kl 



K X K 2 ' 

that in terms of the initial variables reads as 



xJ (x\ — x 2 )x 2 

If we repeat the last procedure but under the assumption that the integral does not depend 
on x 2 we obtain the following first-integral: 

C 2 =[x 3 --)x 2 + — ^2—. 

V xJ (X\ — X 3 )X Z 

It is a long but straightforward computation to check out that both are integrals of M\, M 2 
and M 3 . We can obtain now the general solution x of ([I]) x in terms of x±, x 2 , x 3 , C±, C 2 as 



x _ (Cijxi - x 2 ) - C 2 (xi - x 3 )) 2 + k(x 2 - x 3 ) 2 ^ 
Y {C 2 - Ci)(x 2 - x 3 )(x 2 - xi)(xi - x 3 ) 

where C\ and C 2 are constants such that once x\(t), x 2 (t) and x 3 (t) have been fixed they 
make x(0) given by (|TS|) to be real. 

Thus we have obtained a new mixed superposition rule which allows us to express the general 
solution of the Pinney equation in terms of three solutions of Riccati equations, and of course 
two constants which determine each particular solution. 



5 Conclusions and outlook 

In this paper we have reviewed the theory developed in [6] in order to deal with a certain kind of 
systems second-order differential equations that can be studied as Lie systems, in particular with 
the Milne-Pinney equation and the Ermakov system. We have not only recovered in this way 
some known results about these differential equations, i.e. integrals of motion or superposition 
rules, but we have also been able to show new ones as a superposition rule for three solutions of 
a Riccati equation to give us the general solution of the Milne-Pinney equation. 

New applications of this formalism to relate a Milne-Pinney equation with a frequency uj(t) 
with a ^-dependent harmonic oscillator with different frequency u/(i) related with uj(t) by means 
of certain relations will be developed in forthcoming papers. In this way a new time-dependent 
superposition rule is obtained and by means of it we can express the solution of a Milne-Pinney 
equation with frequency co(t) in terms of solutions of a time-dependent harmonic oscillator with 
frequency u'(t). 

Another remarkable point to be studied is the extension of this formalism for studying some 
generalised Ermakov systems as those appearing for instance in [38] and |39j. 
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